General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



JSC- 09836 


NASA TECHNICAL MEMORANDUM 


NASA TM X- 58166 
August 1975 



A RADIATIVE TRANSFER MODEL FOR MICROWAVE 
EMISSIONS FROM BARE AGRICULTURAL SOILS 


(NA3A-TH-X-58166) A FADIATIVE TRANSFER 
MODEL FOR MICROWAVE EMISSIONS FROM BABE 
AGRICULTURAL SOILS (NASA) 29 p HC $3.75 

CSCL 08M 

G3/4 3 


N76-10552 

Unclas 

03020 




NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
LYNDON B. JOHNSON SPACE CENTER 

t 

HOUSTON, TEXAS 



I. Mftxirt No i Govwnmont Ac cu ion No 

TM X- 581 66 


4 Till* md SubfitN 

A RADIATIVE TRANSFER MODEL FOR MICROWAVE 
EMISSIONS FROM BARE AGRICULTURAL SOILS 


3. Racipionri Catalan No 


ft RapoM Data 


6 Rarformini Orgafiitation Coda 


7 Aullior(il 

William J. Burke, JSC and NRC, and Jack F. Paris, 
Lockheed Electronics Company and University of Houston 


9 Rarforminp Organiiation Nama and Addraaa 

Lyndon B. Johnson Space Center 
Houston, Texas 77058 


13. Soontoring Agancy Nama and Addrau 

National Aeronautics and Space Administration 
Washington, D. C. 20546 


8 RaHofmmg Organitation Raport No 

JSC-00836 


10 Work Unit No 


11 Contract or Grant No 


13 Typa of Raport and Rariod Covaratl 

Technical Memorandum 


14 Sporttoring Agartcy Coda 



16 Abatract 

A radiative transfer model for mlcowave emissions from bare, stratified agricultural soils was 
developed *o assist In the analysis of data gathered In the Joint Soil Moisture Experiment. The 
predictions of the model were compared with pre!. ninary X band (2. 8 cm) microwave and 
ground based observations. Measured brightness temperatures at vertical (T^) and horizontal 

(T„) polarizations can be used to estimate the moisture content of the top centimeter of soil with 

1 1 percent accuracy. It is also shown that the Stokes parameters [T^ + ^"^V " ^H^ 

can be used to distinguish between moisture and surface roughness effects. 


17. Kfy Words (Suggstttd by Author(t) I 

Microwave remote sensing 


Multifrequency microwave 
■ Passive microwave observations radiometer 

* Radiometric microwave observations ’Radiative transfer of 

'Soil moisture microwaves 

‘Joint soil moisture experiment 

* Passive mi crowave Imaging; system 

19. Sacurity Oitsif. (of this rtporti 20 Sacurity Citssif (of this pagtl 

Unclassified Unclassified 


18 Distribution Stttamam 


STAR Subject Category: 
43 (Earth Resources) 


21. No. of Ptgts 

27 


22 . 

$3.75 


’For sal* b>> tha Natiorral Tachnical Information Sarvioa, Springfiald, Virginia 23161 


NASA — JSC 


ItC Fi.V. 1474 (*•« lul 74) 




















CONTENTS 


Section Page 

SUMMABY 1 

INTRODUCTION 1 

THE SOIL MOISTURE EXPERIMENT 2 

MICROWAVE EMISSION FROM SOILS 2 

A COMPARISON WITH PRELIMINARY PASSIVE MICROWAVE IMAGING SYSTEM 

(PMIS) DATA T 

AN ATTEMI^T AT INVERSION 10 

APPLICATION OF THE MODEL TO MFMR DATA IT 

CONCLUSION 20 

APPENDIX — PROGRAM FOR BRIGHTNESS TEMPERATURES 22 




iii 


TABLES 


TalDl° Page 

I PMIS DATA FROM LINES U AND 1 ON THE FIRST FLIGHT 

APRIL 5, 197^ 8 

Ila AVERAGE FIELD MOISTURE AND TEMPERATURE PROFILES ......... 9 

IIL PREDICTED TEMPERATURES BY THE MODEL FOR TOP AND BOTTOM OF 

FURROWS . 9 


III AVERAGE VALUES: G T,, - T„ AND (T,, T„)/2 FOR FOUR FLIGHT 

V TuTtnr* V il V n 


FIGURES 


Figure Page 

1 Stratified medium used in development of emission model U 


2 Emission model predictions for Field 260A. The ascending curve 
is for vertical and the descending curve for horizontal polar- 


izations. Straight lines marked T^ and T^j are PMIS 

observations 11 

3 Model calculations for vertical and horizontal polarizations at 

0^ = 0°, 30° and 50° as a function of the percent moisture 

in the top layer of soil . . . 13 

4 Difference between predicted vertical and horizontal polariza- 

tion brightness temperature for 0^ = 20°, 30°, 4o° and 50° 


as a function of the percent moisture in the top layer .... l4 

5 The average of the predicted vertical and horizontal polariza- 
tion brightness temperatures for 0^ = 30° and 50° as a fiuiction 


of the percent moistvire in the top layer 15 

6 The difference between predicted brightness temperatures at ver- 
tical and horizontal polarization as a function of their 
average for 8o = 30° and 50°. The straight lines are lines 
of constant moivsture in the top layer ............. l6 

T The predicted brightness temperat\ire contributions from the 
first, second, third and fourth layers at 6o = 0° as a 
function of the percent moist^are in the top layer l8 


// _ 


iv 


Figure 

0 Predicted brightness temperature contributions from the first 

second and third layers as a function of the percent moisture 
In the second layer. The lock angle Is 0°. The percent mois- 
ture In the top layers are also Indicated 


A RADIATIVE TRANSFER MODEL FOR MICROWAVE 


EMISSIONS FROM HARE AGRICULTURAL SOILS 

William J. Burke,* and Jack F. Paris^ 
Lyndon B. Johnson Space Center 


SUMMARY 


A radiative transfer model for the emission of microwave radiation from 
agricultural fields was developel. The predictions of the model are found to 
be in good agreement with prelimiiary data from the Phoenix Passive Microwave 
Imaging System observations. A sun angle effect observed in the L band data 
mokes a simple comparison of the model with observations impossible at this 
time. An inversion -of the model was capable of predicting the moisture content 
of the top centimeter of soil with an accuracy of approximately 1 percent. The 
model has been used to distinguish the microwave signatures of smooth and 
rough surfaces . 


INTRODUCTION 


The purpose of this memorandum is to present a theoretical model adi- 

ative emissions from bare agricultural fields. The model utilizes the “ ‘ative 
transfer equation to calculate emissions from stratified soil in terms of bright- 
ness temperatures at vertical and horizontal polarizations as a function of 
look angle and moisture and temperature profiles of the soil. The model was 
developed to assist in the analysis of data from an aircraft experiment con- 
ducted during April, 197^. Data from the experiment have not been completely 
reduced. A small fraction (~1 percent) of the data have been reduced using 
by hand methods and are presented in a preliminary form. The data are used to 
verify the model and point out the capabilities and limitations of the model 
for analyzing radiometric data. 

The first section gives an overview of the experiment EUid is followed by 
the development of the model and a comparison of its predictions with prelimi- 
nary X bemd data. On the basis of this comparison aii inversion is presented. 

A few words of caution are given concerning applications of the model to L band 
observations. A program listing of the model calculations is in the appendix. 


*NASA Research Center. 

^Lockheed Electronics Company and University of Houston. 



THE POT I, MOISTIIRE EXPERIMENT 


A feoalbllity test for detecting soil moisture using microwave remote sens- 
ing techniques was conducted under NAPA auspices near T’hoenix, Ailrona April 5 
and 6, 197*<. Contributing Investigators come from the Johnson Space Center, 
Goddard Space Flight Center, Agriculture Research Service of USPA, Environmental 
Research Institute of Michigan, University of Arkansas, University of Kansas, 
and Texas A&M University. X and L band radiometers were fh^wn over four flight 
lines. Lines 1 and ? are 3? kilometers (20 miles) long and alined porth-south 
and lines 3 and L are 6L kilometers ( Lo miles) long, east-west alined. Three 
flights were made over the flight lines during the early afternoon of April 5 
and a single flight at dawn April 6. At the time of the flights soil arid tem- 
perature groiuid ti'uth sanpleu were token from 9^ l6 hectares (forty acre) fields. 
Moisture samples were token from the tops and bottoms of furrows a* depths of 
0 to 1, 1 to 2, 2 to 5, 5 to 9 and 9 to I'' cm at four points in each field. 
Temperature samples from the center of each layer were taken at one site per 
field. Soil samples were towen to the Agriculture Research Service La*^ ratory 
in Chickasha, Okalahoma for hygroscopic (ref, l) ojid texture analysis. To 
measure the dielectric coefficients at X and L band frequencies sarai'les of 
representative soil types were taken to Texas A4M University. Tc date, the 
dielectric coefficient analyses at X band have been completed but the L band 
have not. 

The X and L band radiometers carrie>l aboard the NASA P3A aircraft are the 
Passive Microwave Imaging System (PMIS) and Multi frequency Microwave Radiometer 
(MFMH). Tlie PMIS (ref. 2) is a 10.69 GHz Imaging radiometer that scans in I 4 I 4 
steps along an azimuthal arc of *3l*.5° at a constant look angle of 149.5°. The 
antenna is a dual polarized cross slot array electronically stepped for scan- 
ning. Beef 'se of })rogrnmming difficulties, PMIS data are not ccwipletely re- 
duced. At the time of the April mission only the L band (I.L 2 GHz) component 
of W’MR was Installed in tlie aircraft (ref, 3). ITie antenna is a flat plate 
dipole array whose look angle and polarization must be manually set. On the 
first afternoon and on the dawn flights the antenna was set at nadir look. The 
look angle was L0° for the second and third afternoon flights. Vertical and 
horizontal polarized brightness temperatures were rieasured on the respective 
flights. The M^MR data have been reduced. 


MICROWAVE FWISSION FROM SOILS 


This section discusses the problem of calculating the intensity of micro- 
wave emission measured by an antenna above an agricultural surface. Developing 
the radiative transfer model the following simplifying assumptions were made: 

1. The radiation is incoherent 

2. There is no attenuation or emission between the surface and antenna 

3 . Tlie sky brightness is isotropic and has a value of 30° 
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U. Moisture tind temperature are functions of depth only 

5. Dielectric properties are constant across any given layer of soil 

6. The surface of the soil is smooth 


The extent these aasump^ ions brea,K dowi* is discussed where the predictions of 
the model are compared with PMIS and MFMR data. 


A cross section of a stratified soil is shown on figure 1. Layers have 

thickmsses AZ , which are not necessarily the same for all layers. The J th 
J 

layer is banded on the top by the J th surface and by the J th surface on 
the bottom. Within this layer the dispersion relation for electromagnetic wave 


propagation is k 


J 


(-) 

c 




The frequency is w in radians/sec, c the 


velocity of light, u the magnetic permeability (assumed equal to one) and 


'J 


♦ ic is the complex dielectric coefficient. Written as 


RJ 


IJ 


^ ^hen 

J C J J 
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Snell's law shows that the component of the wave vector parallel to the surface 
is a conserved quantity. Thus 


2 2 
6 , + 6 , 
xj yj 


constant “ Sin 0 o 


and 


a , = a , 

xj yJ 


0^ is the angle the ray emerges fr.m the soil with respect ti the surface 
normal. It ranges in value between 0 and 90®. Equation (l) reduces to 


and 



c - Sin ^0o 
Rj 
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IJ 


( 2 ) 
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Solvliig 



At the boundary between the J and J-1 layers radiation is partially reflected 
and transmitted. The fractions of the incident - iectric field with horizontal 
and vertical polari zatioi.o reflected back into the J th layer are given by 
the Fresnel coefficients. 
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CJ-l k^, - CJ k^^_^ 

^ CJ-1 k^j . CJ k^j.^ 

where k , ■ S , ia . Thus 6 a , p and p„ all depend on the complex 

ZJ ZJ ZJ Z Z il V 

dielectric coefficient and angle 0^ that the ray emerges from the soil. 

An attempt to construct a radiative transfer equation that describes radi- 
ation emitted from a stratified soil is now made. Within the first layer the 
radiative transfer equation is 


d I 
d 2 


I + J 

(i) 1 U) 


(5) 


Y is usually written as a product of the density and monochromatic mass absorp- 
tion coefficient. By writing the Poynting theorem in an appropriate form it can 

be shown that y, = ?wa ,(0 )/C. 

i zi o 

I is the intensity of rad.^ation ct frequency w. J is the Planck 

Ul U) 

emission function. In the microwive frequency range Planck's emission law re- 
duces to the Rayleigh-Jeans approximation, where J is proportional to the 

U) 

temperature of the rntdixm T. Adopting a similar scaling rule for I^, and 
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effective temperature T can be defined which is directly proportional to I . 

p (1) 

The buoscript w is suppressed and refers to the intensity in a narrow 

riuige near u. The subscript p refers to a polarization state. Since J 

U) 

is isotropic and independent of polarization no designation is necessaiy. The 
radiative transfer equation in the first layer may be written 


d 




T, 


( 6 ) 


This equation can be integrated from a point Just below the surface to a point 
Just above the interface between the first and second layers. Because the 
dielectric properties are assumed to be constant across the layer 


T (l‘ 
P 




(2*)e 


-Y Az 
'l 1 


( 7 ) 


The argxunent Ni implies that the measuroment is made above (•►) or below (-) 
the N th interface. The fi”st term on the right hand side of jquation 7 ac- 
counts for radiation emitted within the first layer and comes directly to the 
surface. The second term describes upwelling radiation at the bottom of the 
second layer. This in turn has two components: first, radiation emitted In the 

first layer and reflected at ^he interface between the first and second layers; 
and second, radiation transmitted from lower layers. 


T 

P 


(J*) = Rpj Tj 



-Y^Az 




+ T ( 2" ) 
P 



( 0 ) 


B „ is the absolute value squared of the Fresnel coefficient for p polariza- 
P2 

tlon (equation Ji). The radiation field Just above the surface is the value 
Just below multiplied by the transmittance 


T d"^) = (l-R l'^ T (T) (9) 

P V P / P 

■ ('-v) ! 
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The radiative transfer equation can be integrated again to calculate T (? ). 
Repeating the procedure for N layers gives ^ 


T (!♦ 
P 


,e ) • 7 T. \l-e ^ ° / IfR (e ) e 

O i Pti+l o 

i-1 

i 
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\i (10) 


n (l-R (6 )) <• 

P.J o , p 
J-1 




The brightness temperature measured by the antenm. is the sum of the reflected 
sky brightness and radiation emitted from the soil 



(Q ) ■ T 
o 


SKY 



(O ) ♦ T 
o p 
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A computer program was written that used as its inputs the moisture and tem- 
perature profiles measured in the Phoenix fields. Dielectric coefficients were 
calculated using the soil type characteristic of each fielci with the polynomial 
representation given from Texas A&M measurements. With this information 
equation 11 could be calculated and results of the model compared with PMIS 
observations . 


A COMPARISON WITH PRELIMIl^RY PASSIVE MICROWAVE IMJUIING SYSTEM (PMIS) DATA 


April 5 and 6, 197^* missions PMIS data have not been reduced to the point 
where average temperatures and standard deviations over fields are readily 
available. However, it is possible to take data from uncorrected printouts 
and estimate average temperatures by considering only data taken near the 
center of a field. The preliminary nature of average values cannot be over- 
stressed. They do serve the useful purpose of helping to check the model's 
validity. The PMIS data are taken from lines U and 1 on the first flight of 
April 5 flight, and are given in Table I. 

A program was developed, which tfJtes the observed moisture and temperature 
profiles as well as soil types and calculates the horizontal and vertical 
rightness for the top and bottom of furrows in ten degree increments fiom 0 
to 90°. A quick preliminary check of the specular model's applicability is to 
compare whether the predicted temperature at a look angle of 0° falls between 
the PMIS Tv and Th observations. No field failed to meet the minimal criterion. 
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TABLE I.- PMIS DATA FROM LINES U AND 1 ON THE FIRST 


FLIGHT APRIL 5, 197l« 


Field 

Tv 

Th 

PM '1) 

237 

287 

273 

1.6 

2 U 2 

288 

275 

1.6 

2^3 

286 

21U 

1.7 

25 I 4 

289 

27h 

3.3 

255 

291 

27h 

U.l 

257 

.'86 

269 

2.1* 

260 A 

283 

267 

I4.I* 

261 

?8t 

272 

3.3 

26U 

286 

271 * 

2.8 

299 

281 

266 

2.5 

300 

285 

266 

2.3 

317 

288 

251 

1.7 

33 U 

28!4 

272 

1.1* 

102 

277 

266 


lOU 

283 

271 


106A 

283 

261 

8.8 

126 

285 

27h 



Field 26 OA was selected to illustrate the capabilities of the model. 

Ground obsei-vations were simultaneous with the overpass. Avera^^e moisture and 

temperature profiles for the field are given in Table Ila. There is a wide 

difference between the moisture in the 0 to 1 cm layers of the top and bottom 

of furrows. Temperatures predicted by the model for the top and bottom of 

furrows are given in Table Ilb. The temperature 272® K at 0® look angle from 

tne top of the furrow lies between the PMIS T„ = 283® K and T„ ® 267® K 

V H 

observations; the 225® K temperature of the furrows bottom does not. It is 
concluded that most of the radiation reaching the PMIS antenna comes from the 
top of furrown. Furrows in field 260A had a height of ~27 cm at a separation 
of 1 meter and were alined perpendicular to the flight line. Thus, significant 
shadowing of the bottom of the furrow 5s not a complete surprise. 
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"ABLE Ila.- AVERAGE FIELD MOISTURE AND TEMPERATURE PROFILES 


Field 26 OA 


Layer 


Moisture 

Temperature (®K) 

u.u 

305.5 

10.6 

303.7 

19.7 

301. L 

21.6 

297.3 

23.7 

290.3 

12.8 

30 I 4 .I 4 

17.5 

302.0 

21.1 

299.2 

23.6 

295. 

23.0 

289.3 


271.9 

?72.h 


275.^ 


280 . L 

287.3 

295.6 

302.6 

299.3 
253.5 

30.0 


71.9 

71.2 

68.2 

62.6 


253.2 

238.3 
21L.7 
177.5 

119.3 

30.0 


Bottom 

'IV 

Th 

225.2 

225.2 

225.7 

22L.9 

230.2 

220.5 

238.0 

212.6 

2 U 9.5 

200.5 

265.0 

181.2 

283.9 

159.3 

300.8 

127.1 

285.3 

8 L .7 

30.0 

30.0 
























r>’edicte<<. temperatures from the top of furrows as a function of look angle 
are plotted or Figure 2. The straight lines marked T,, and T„ are PMIS observe- 

tions. These curves intersect near a look angle of 30°. An inspection of data 
from the other fields show that with the exception of fie^d 317* the intersec- 
tion of predicted and observed temperatures occurs near 30°. Field 317 was a 
flat, smooth field. Here PMIS observations intersect the prediction curves near 
b*^0°. The conclusion is that surface roughness changes the effective look 
angle from 50 to 30°. Such an effect was anticipated by predictions of a rough 
surface emission model develop'- 1 by Ulaby et al. (ref. U). Their geometrical 
optics model peedlcts that ns roughness increases the vertical emlssivlty is 
depressed and horizontal eralssivity enhanced relative to the emissivity of r. 
specular surface. Thus from a specular point of view, observations made at a 
look angle of 50° appear, to be made at a lesser angle. The apparen*. univer- 
sality of the 30° observation suggests the roughness model can be Inverted to 
determine the size of eleds responsible for the X band emissions from Phoenix 
fields . 

'TVo sets of calculations were made using varied layer patterns. For the 
first calculation temperature and moisture of the first two layers were 
averaged and combined. For the second calculation, all layers were set at a 
thickness of 0.5 cm. The temperature ana moisture in each layer was set as the 
average of a linear interpolation across the layer; neither satisfied the zero 
look angle criterion as well as the oi-dinary pr'iile. 


AN ATTEMPT AT INVERSION 


If a model is to be useful in remote sensing it must be capable of Inver- 
sioji. Wo must be able to go from aircraft observations to an estimate of the 
soil moisture. To effect such an Inversion, the brightness temperatures of a 
set of pseudo- fields have been calculated. The temperature profile from top 
to bottom of these fields was set at T^^ = 303, T^ = 301, T^ = ?99* T^ = 297, 

T, = 295. The moisture profile was PM(l) = M, PM(2) *= M + 2, PM(3) = M + J, 

PM(l4) - M + II 4 , PM(5) = M 16 . Tlie dielectric coefficient is that of sandy 
clay. M was given values of 1 to 25 but the percent moisture in any given 
layer was not allowed to exceed 30 percent. This moisture profile is referred 
to as the ordinary profile. Other moisture profiles with sharper gradients 
were used. Unless otherwise specified we will be referring to model results 
using the ordinary profiles. 
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Figure 2- Emission model predictions for Field 260A. The ascending curve is for vertical 
Ind the descending curve for norizontal polarizations. Straight lines marKed Tv and Th 
are PMIS observations. 




Values of calculated brightness temperatures as a function of the percent 
moisture in the first layer PM(l) are plotted on Figure 3 for look angles of 
0, 30, and 50®. At the 0® look angle and PM(l) > 10 percent the average slope 
is -3.5® K/percent moisture. The slope tends to decrease at higher values 
of PM(l). For this reason it is thought that the model's predictions are con- 
sistent with the observations of a -3.3® K/percent moisture slope reported by 
Schmugge (ref. l) for moistures up to 35 percent. Also the plot on Figure 3 can 
be used to attempt an inversion. Unless the presence of a flat field is indi- 
cated by the large separation between and Tj^, it is assumed that PMIS obser- 

vations come from an effective look angle of 30®. The estimates gained by 
this method are found to vary by no more than 1 percent from the average values 
of PM(1) found in ground truth results. 

The predicted difference between T^ and T^ as a function of PM(l) at look 

angles of 20, 30, Uo and 50® are plotted on Figure 1*. A scatter plot of FWIS 
observeu values of T^ - T^^ as a function of observed IW(l) is superimposed cn 

the calculations. Average values ">f fields 317 A, B, C, D ore found near the 
0^ ■ 50® curve while the lest of the data points lie close to the 0^ ■» 30® 

curve . 

Figure 5 is a plot of the average value of the vertical and horizontal 
brightness temperatures as a function of PM(l). The curves are for 0^ * 30® 

and 50°. We note that 1/2 (T,, ♦ T„) is only slightly dependent on the surface 

V n 

roughness characteristics. A scatter plot of PMIS observations appears less 

correlated with the 1/2 (T„ T„) predictions than with the (T„ - T„) predic- 

. , V H V H 

tions . 

A plot of (T,, - T„) as a funciton of 1/2 (T,, T„) for 0 “ ?0 and 50° is 

V H V H O 

given on Figure 6. The nearly vertical curves are lines of constant moisture. 

A scftter plot of PMIS observations is also displayed. Data points taken fro.u 
fields 313 and 296 appear on the graph. Aerial photographs show the fields were 
being irrigated d’oring the overpass. No ground truth for these fields la avail- 
able. Th^ anomalous positions of the points on the graph suggest that the beam 
is made up of components from smooth water flowing in furrows and relatively 
dry tops of furrows. Signatures from a small pond and large metal building on 
line are plotted. The pond is cool and flat but the metal building is cold 
with no appreciable polarization. 

It should be pointed out that if 1/2 (T„ + T„) and (T„ - T„) should turn 

V n V n 

out to be significant paramete.’s, it would not be a total sui'prise. Paris 
(ref. 5) has shown that these quantities are the non-zero components o<' the 
Stokes vector. Effectively the model sums up the contributions of the Scokes 
vector for each layer. 



PM(1) 


Figtire 3.- Model calculations for vertical and horl ?.onta!! polarizations at 
Go = 0®, 30° and 50° as a f\inction of the percent moisture in the top 
layer of soil. 







PM(1) 


Figiire 5-- The average of the predicted vertical and horizontal polarization 
brightness temperatures for 6o = 30® and 50° as a function of the percent 
moisture in the top layer. 
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Fieure 6 - The difference between predicted brightness temperatures at vertical and horizontal polari- 
Ltion'as a function of their average for Oo = 30° and 50°. The straight lines are lines of con- 
stant moisture in the top layer. 





The model provides a means of knowing from what depth the radiation reach- 
ing the antenna comes. Figure 7 gives the contributions fnjm the first, second, 
third, and foiirth layers aa a function of PM(l) for the ordinary profile. The 
contribution from the first layer is greatest for PM(l) ■ 11 percent. Although 
the radiation generated in this layer increases with increasing moisture, the 
fraction allowed to escape decreases. Output calculations were made with dif- 
ferent moisture gradients. Contributions from the first three layers as a func- 
tion of PM(2) are plotted on Figure 8. Beside each curve is the percent mois- 
ture in the first layer. Although the contributions from the first and second 
layers Incre'se with increasing PM(2), the output from the third layer mono- 
tonically dec.'»ases. The sum of the contributions, except in the case of very 
dry crusts, » .lows little information about the actual gradient that gets 
through. In fact it is impossible to distinguish between cases in which 
PM(1) * 1 percent, l’M{2) * 20 percent from PM(l) ■ I 4 percent, PM(2) ■ 6 percent 
even though the average moisture over the top two centimeters is quite different. 
For PM(1) > 10 percent information abc ut lower layers is completely lost. Per- 
haps, the ambiguity noted by Schmugge (ref. l) for average moisture of < 10 per- 
cent can be partially explained in terns of the model's predictions. 


APPLICATION OF THE MODEL TO MFMB DATA 


Because the dielectric coefficients at L band frequencies of all of the 

soil types of the Phoenix area ore not available, it is impossible to compare 

model prediction directly with MFMR data. However, an indirect attempt has been 

made using available T„ and T„ data for the April 5 flight. This attempt was in 

V H 

the form of a scatter plot of T„ - T„ as a function of soil moisture. Prelimi- 

V H 

nary results we»'e not encouraging. Part of the difficulty arises from the fact 
that the aircraft flew north on line 1 and south on line 2 while making 
Ty measurements. The opposite directions were flown while making T^^ observa- 
tions. Thus on line 1 we observed the cold side of ridges measuring T^ and hot 
side measuring T^^, and vice versa on line 2. The effect is to diminisii (T^ - T^^) 

observation on line 1 and enhance them on line 2. Sun angle effects on the 

parameter (T„ + T„)/2 should cancel. 

V H 

As an example of this sun-angle effect consider observations made in 

Fields 11 and 13. The soil of both fields is sandy clay loam. The moisture 

profiles were similar and state of cultivation of the fields was identical. 

The temperatures observed in field 11 (line 2) were T,, * 292.1 and T„ = 273.3, 

V n 

while for field 13 (line l) T„ = 269.6 and T„ = 276.1. A similar systematic 

V 1a 

dimunition and enhancement of (T - T ) is found in all of the lines 1 and 2 ob- 

V n 

servations. Evidence for a sun angle effect is also obvious in the lines 3 and 
I 4 data but to a lesser extent due to the fact that the midday sun was shining 
almost directly c it of the south. The average values of T^ - T^^ and (T^ Tjj)/2 

for the four flight lines (Table III) show an enhemcement of the polarization 
difference along lines 2 and U and a dimunition along lines 1 and 3. Within 
a standard deviation, the intensity is the same along the four flight lines. 
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Figure 7.- The predicted brightness temp>-rnture contributions frn t » » 

second, third and fourth layers ~ 0® as a functloi. t f •’ •)> r -eni 

moisture in the top layer 
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Figure 8.- Predicted brightness temperature contributions from the first sec 
ond and third layers as a function of the percent moisture in the second 
layer. The look angle is 0°. The r'>rcent moisture in the top layers are 
also indicated. 



TABLE III.- Ay^ERAOE VALUES: 0 Ty - '1'^^ AND (’’.y ♦ 

FOUR FLIt’HT f INcS 


Line 

T - T 
W H 

(T., ♦ '1 'fJ 
y t! 

1 

9.^ t 2.^ 

•’03.8 ' . . 

2 

1W.6 ‘ L.4 

279.'* * 9.y 

3 

11.2 ? 7 

,^83.8 * 3.5 

U 

,0.7 1 5.2 

281.2 ♦ S.b 


Another pertubation comes l>om the scales of roughne.''^. It is the 
relatively large facets that con’^ributc most t& X btind -mi . on Pending 
further analysis of the roughness scales, it seom.': (ju’le i eascnuole to expect 
that large scale roughness at X hand may turn out to be sinal’ an 1 mi'ldle stale 
at L band. This is not to say that tfie mo iei ir inap] 1 1 tabi c, rather that it 
has not been applied. Also the application in the analysis of data tiay be far 
more subtle than the X band application. 


CONCLUSION 


The radiative transfer e<juation to the problem of micr'jwave einibsioiis from 
layered, moist soils were applied. The predictions of the mole' «^re conp’reo 
with a small set of X band observations and it was shown 'ha! • 

1. Except for the case of flat surfaces the effect of s^ifac» jourhuers 
is to change the effective look angle from ^0® to to'' 

?. The predicted slope for nadir look is consistent with | reviiin ^ . — 
ported observations (Schmugge et al, 197*t) 

3. The model can be Inverted to give the percent moisture in tfi« - p 
layer 

'Hie third conclusion does not preclude the possibility of gaining i nfoj': i < 1 Ion 
about moisture in lower layers. We doubt that it can be known with the n '-u acy 
that we have of the top layer's moisture content. 

One final result which deserves some stress concerns the role of th< 

parameters 1/2(T,, + T„) and (T,, - T ). We have seen that these arc tie ''om- 
V H V H 

ponents of the siunmed, radiatively trans'’erred r^tokes rect.c*':^. rhe d 'St»'ibut ic.. 
of points on Figure 8 suggests '.ha*, trice • a’''urerers can bf- urei *o make 



maxlimiffl use of polar! zatior information ncutained in the P4IS data by way of 
field classification. A cursory look of imagery suggests that these 

parameters can be exploited to classify various agricultural scenes. 

The paper has not presented our model cs new. Rather using well estab- 
lished principles, it has suggested heuristic, i^ethods by which information con- 
tained in microwave data may be exploited to ?i h*vnce our understanding of this 
remote sensing technique. 
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APPENDIX 


PROGRAM FOR BRlGirr :E8S 'nWP)-JDVTtJRt'S 


OlMLNklOM 

OINCNSION 


OlbRCtllOit TmIjIOI, TM^j iOl ,T¥l , «0| |TVi« lOl 


IPSt ( 7 *a> | 6 | 7 • 2 ^ • 


I «<IO>7.2*ifcllO*7i2i«llloi,VllC*/.2) 

0 I M| NS I ON «l ICI 7 |2I .SUN I t0.7.2»,*RH<IOi7i2tiTNvilUi7i2l,T,H(IO»7L2 
1 I • IgV( I 0 I f t Z i iRh^ 10.7 I Z t IRVl I0«7.^) 

COMfLl> crsl7,2). ^»VI|0|7.2I, *hSh * i 0 . 7 , ^ i • NhuVII0i7i<) 

InTcoi" sano 


1 • lfc» 0 <* I *t*/ 

OA^t OZ* I * 'OZ'^* lOZ*^* lOZ**** »02*S* »OZ*A* .^Z*7l/0»OiJ»(‘.J»O»^«0i 

I H.0iA*Oit*lA^ 

OaTa P I ' J* I •• I S72*Sn^ 

HtAO Ib.lbOOlBANO 
I 000 7 0**NA T i ( S > 

c t»'7tP ritto iUIL TYPt, NoiSTjBt ayU T A 7uK£ PnOf.tt^ 

NtA0lfcij7MTCT«l»,TcB«l»«|»l,S» 

J7 YONNATIlOfStl) 

1 "tAOIS.ZI 7N ,a5T , H P*"** « • . J J , I • I ,S » , J* I , Z • 

2 FuhnAT I AS I t S • I 0F& . n 


C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


7 N ■ F I NuNbCR 

aST» soil TTRt 

Tnt UUA(I.T|TT PNb |S ThI AVtAAftt "OlStUNt fOl/Nj IN 

•i. AnCNANO • S SMO^NW T<*u7^ SUQK 

Pt*»CtNT NolSTUPt IN |TM LAFt« 
ni.J>« TtNPtBATURt NEASUHCO |N ITH LA^Sk 


J • I SANPut 7*0H TqP of FUiYNO* 

J • 2 SaNPLI f"ON SOTTOM of FuPRO* 

KS7 • I sandy cI^ay Loan 

KS 7 • 2 CLAY COAH 

AST • J S*nOt cl-*T 

K^r • H CLAY 

IF BANU is I AE ABC USC|NS L BAND 2 IS F QN X ^aNu 
initialize NEEPiO NATB|CES 


Ou 1*3 J* I I 2 

©V l »3 ^■l ♦7 

Oo I B) A* I I I 0 
Pn I I • U> * 0*0 
HI Til . Jl"0»0' 

ThI IK)«o>0 
T h2 I K I •0« 0 
.. . . txUJUPXUO 


w pace J 3 


T v2 I K I •0* u 

cpsNi : iji«i«o 
cpsi 1 1 iji«o«o 

tPSII.jl* C«Pl* ItpSRI I I,, .r.r'SI I I I J» » 

SuHIK • I • J> • 0*0 

. ... 2 A _»_-i.t 0 -- - 

TNnIK • I • jI ■ 1 iQ 
IPS ThyIa » I • JI • I to 

BMOHlK.I.JIalOtQiOtol 
BnOVUtl • Ji*<OYOiOtc’ 

BhIK.I.jI • OtO 

-? . 0 x 0 

ThN|!(,|,J| ■ Q.o 

|B3 TbvU, I , ji • 0*0 

C bCC*J&C TH£ soil Tenper^TURES BCRC Be*0 in FRoN BOTTON TO fOP 
C|T Is necessary to InveR^ the TCNPERATu*IE vector# TNC CORRE|TbO 
C vector is SIVEN the NAhC tcct or tccb# 

JiiT IJJ.-1CT1S1 

T C C T I 2 I • T t T I •• I 


TtCT«3»*TCTU| 

TjcT «‘il»TcT»2» 

TccT «H»TcT« J 

'ecu* * 

Uc*‘i'*Tco<*'» 

T t c B ' • T c 9 ‘ » 

UCB< & > •TC9< II 
Ob *•; J*ti2 
h 2 T I 1 . JI ajC'O 
Do -ii I *2 • 6 

ni.li ■nitll *TcCTll-ll * 273.0 
i3 Iili2l • TUtJI ♦TCC«‘l*ll ♦ 273*0 

Do *« * ja 1 • 2 

Uw *1* la2l* 

•1* PhllijIaPhSlI'llJ) 

00 *<a j> I I 3 

7ijlapn(6ij) 

T|7,J|aTI«iJl 
60 TO l2&3(2&llia«Nj 
3 i| 6o 70 « i»sr 

3 Ou 73 jb|| 2 
Do 73 l»2,7 

a* If tPHi I , jl’S.O' I 3.23*23 
|3 tPSm I • J>"2.i»*0.2‘473apMt I , jl 
tPill|.jl*0*‘*2 
60 TO 73 

23 tl'if'* I •j>*i*'^31-a.2‘,’apMl |.j»*0*0*VaPMl |•Jl•‘‘2-0•UUl•^'Ml lij»a*3 
fcPSlI |.j»"2*2oT- 0«6S’*^M* 0*6aPMl lijia'A-O’OOl'I'Ml Iijlaa3 

73 tPilliJl* CMPl* ItPiTI I • J» *tPS| I I I jI I 
60 TO 100 
i Ou 7H jai ,3 
00 7*« Ia3.7 

09 I F ( 7*'^ I I • J I '« • 0 I 1 0 , 3<* • 2*« 

I tP iN I 1 . j I •*: •0*0. I I H3SaPM I 1 , j I 
tPil I I . Jl«0* 7b 
6o TO 70 

<** tPb«< Iij>"-3»i 73*|. iab»PM( l.jl-0*C22*P«l l.j)»a2*0.000JapH| |.j|a«3 

tr'ii' i*^*“OiO30*r*’U*j''* 2*o.ooo!»at«ij .ji __ 

7** tPill.Jl" C«*'L» ItPStTI 1 . jl < 1 . jl I 

6o TO IOC 
5 00 7b J»|.2 

00 7b 1-3.7 

«0 |F(PFil|.ol-9)lbi2t>»2b 

LrSmi.0l*2*b*Q» 32b *>»H(l.j) 

tPb|(l.dlaO*b*0*0»a^"‘l*jl 
60 TO 7b 

2b tPSmi.J»a-4,iaUl,7bBaPM(i,j|.o.Ob2aPMl|,dlaa3*O.OOlae"»i*J)**-* 

t^SI I I .jI"“*. 33I-1 .2*0a^hl I .d'-0.0**9aPMl | ••3a0*0009a^n« 1 ,jla*J 

7b tPbll.dia CHPlX ItpbRI I • ji .EPS) I I . J> I 

_ _-40 15 105 

2bO Oo V I oa| .3 
00 »l |a2|7 

EPS9I I .dl«2.7|2-0.9o^**^" « I » J* - 0 . I 3 I 3 • P" (l*j|a*2 
1 aO*007*l’»PP I I t j I a*3-0» 000 I 1 7aP(nl j I 0 I ••a 
tPSl I I . jl*-0*00’7b*0*2<«9*apM I ] . j I -O • 03b*9 0 PH JI.jl**2 

1 ♦Q. QQlt9 2*PH ll. j)a* 3 -0 . 0QQ02 7 2i • pH 

91 EPSll.J)>CHPLxlEPsRlltJ>lEPSMl.O)l 


1 . 


r\r> r> 


In this Part of The pkog«*K/ the values of ihe real ano iHaiiINa«t 
IPaRTS of The *aVE vector aS AELL aS The extinction cOfc^fltlENT 

2*Rt calculated as a fUNcTjON Or The Angle uF tHEf’tENtti- 

100 u*l»2 

Ou 6H t«l ,7 
OU 6H K« 1 • 10 

OEOREeU* ■ FlOATJK- 1 J • 10 
X(X)"OEGHtE(K)*Pl/iao*D 

TlK,I,j).tPSH(l,J|-SlNlX<»Cll.,2 - — - 

BU»liJ»"SaRTI0#5*Y(K*liJ)«ll-»SBRTl !♦< EPSa l»j}/lyUi|ij»J)**2Ml 
AlKHiJ)*tPSIIiiJ)/|2>0«li(KtIijn 
IF (6AND-1) in«l|t,n2 

111 0 (X , I , J) b0.6*A ( K • 1 , J) 

GO TO A<( 

112 GIK. 1 ,01 .HB*A« X, I I - - - 

6<4 CoVU,l,j) >CHPLX (blK, I , Jl • A(KM iO> 1 

00 261 J-lt2 
Ou 26 1 !•! |6 

Ou 26 1 K«>1H0 

HHOH<X»I*J’*‘c*'V«Xtl*Jl-CYlV«Xi|*l,Jl)/(c'*ViX*l»J>*C»'^U|l+liJ>) 

RhUtltj) • CAOSlRHOHUt I I J» >««2 

RrtOv«X,I,JI-lEPs(liJ»*CrtV(K,l*i,Ji-EPsU<-|,dl*CAVlR»IiJ>l/(EPs(l* 
lJ)«C«VV(X.itl,JUEpS(UI»Jl«C«V(K»l,J)| 

2ot KvUilij) * CaBSIrHuVUi I • Jl )**2 
c calculate The TRAn&hISSIOn factor FQR N An^ v polarizations 

c Calculate the weighting functioni a(Xi1|J» 

Ou 8 1 d-l |2 - 

DO 8| l>2t7 
Ou 8| KbI >10 

DU 81 Hb2 • 1 

TRhIX>I»J) ■ TRHIK, UC - Rrt‘K»H-l*Jl> 

82 iRVixiliwI) • TrVIK, I>Jl*n*0 -RVU,M-1 , Jl 1 

SuMtK»l>jl • SUM<K*|'U* * GU.H-I.JI'dZIH-iJ .. 

81 aU»IiJI “ EXpt-SuM(K» li J1 1 

C calculate the BR 1 GHTNESS ■ temperature 

DO S3 us 1 , 2 
00 u3 KbI ,13 

Tbril K , 1 , J > • Til , J) •rHU , 1 , J» 
o3 ToVlK,l,J»« T( I ,J1 *rV1K, I , J» 

DO 3H J» I >2 
00 8H 1-2,7 
0 0 8 M K » I , 1 0 

TdH(X,|,j)» Teh(K,l,J» ♦ T U »U> •»! X, I , J> ‘ThhI X* I t Jl • < 1 .0-tXPl -G1 A • 

I 1 » J > *02 < I I M * « 1 , 0*RHl X , 1 , J1 BE*P I -G( Xi I , Ul 'Oil I » > > 

TbV(K, 1,J)« TttVlK,l,,J) ♦ tI 1 > J» ‘WIX, I , J) *ThVI X, 1 > J)« < 1 ,JJ-LXP< -61 
I I > U I *02 ( 1 1 H • ( I • 0*R v( K , I , J 1 BEXp < -G1 K , 1 , u) "U^l 1 J) 1 

CONTINUE 

00 7 9 X- I , 10 

1 « I 

lb9 ThIU ( -THU X ) bTbH(X , 1 I I ) 

Trt2U)«TH2(XUTBHiK, I ,21 
TV1<XI-TV11 XUTbV(K, 1,11 
Tv2(X)BTV2(Xl*TeVlK,l,2> 

I ■ 1 

1FU*LE,7) GO TO lb? 

99 continue 

WRlTEI6,33l ) 

331 FURMATllHl,’ F.N XST BaNO Pm f“ 


OKIGINAIi PAGE US 
OP POOR QUAIiCBfl. 
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*NlTU6.Mb0irK»KSTiiiAN0 

H !>0 r umhaT I • 2 (• I 
*X 1 I 6 , I 3 J I 
IJJ KoNMaTIIAi//* 

AnlTCiAiNis^ PHfT 
Hab >0*<HAT‘yFli»Vl 
mKI r^(6it*30> 

• JO FuMnAT<|Hli' Tf,H TiiV 

AH|T£IA,lJ2i 
J J2 ^UHMAT J 1 a I // ) 

ANlTC.(t,eiC> TBHiTBV 
to I 0 f UKM at ( 7 ( I Ot 1 I • ) / » 

AMITEI6, tbBI 

li»8 F0HHAT«|Hii* Thl Tm 2 T^l 

AKiTEItoi IJ**) 

J JS r OHrtAT » I A . // » 

i*k1TU6,71»THi,th2,iV»,TV2 
71 Format (10eIi*s/) 

Go TO I 
L i< U 


• . » 


Tvi • i 


c M p I u * T I 0 ^ ; 
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